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Abstract 

One of the most important problems in astrophysics concerns the nature of the 
dark matter in galactic halos, whose presence is implied mainly by the observed 
flat rotation curves in spiral galaxies. Due to the Pauli exclusion principle it can 
be shown that neutrinos cannot be a major constituent of the halo dark matter. 
As far as cold dark matter is concerned there might be a discrepancy between 
the results of the N-body simulations and the measured rotation curves for dwarf 
galaxies. A fact this, if confirmed, which would exclude cold dark matter as a 
viable candidate for the halo dark matter. 

In the framework of a baryonic scenario the most plausible candidates are brown 
or white dwarfs and cold molecular clouds (mainly of H 2 ). The former can be 
detected with the ongoing microlensing experiments. In fact, the French collabo- 
ration EROS and the American-Australian collaboration MACHO have reported 
the observation of altogether ~ 10 microlensing events by monitoring during sev- 
eral years the brightness of millions of stars in the Large Magellanic Cloud. In 
particular, the MACHO team announced the discovery of 8 microlensing candi- 
dates by analysing their first 2 years of observations. This would imply that the 
halo dark matter fraction in form of MACHOs (Massive Astrophysical Compact 
Halo Objects) is of the order of 45-50% assuming a standard spherical halo model. 
The most accurate way to get information on the mass of the MACHOs is to use 
the method of mass moments, which leads to an average mass of O.27M . 
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1 Introduction 



One of the most important problems in astrophysics concerns the nature of the dark matter 
in galactic halos, whose presence is implied by the observed flat rotation curves in spiral 
galaxies [[l], @], the X-ray diffuse emission in elliptical galaxies as well as by the dynamics 
of galaxy clusters. Primordial nucleosynthesis entails that most of the baryonic matter in 
the Universe is nonluminous, and such an amount of dark matter falls suspiciously close 
to that required by the rotation curves. Surely, the standard model of elementary particle 
forces can hardly be viewed as the ultimate theory and all the attempts in that direction 
invariably call for new particles. Hence, the idea of nonbaryonic dark matter naturally enters 
the realm of cosmolgy and may help in the understanding of the process of galaxy formation 
and clustering of galaxies. 

The problem of dark matter started already with the pioneering work of Oort in 1932 
and Zwicky || in 1933 and its mistery is still not solved. Actually, there are several dark 
matter problems on different scales ranging from the solar neighbourhood, galactic halos, 
cluster of galaxies to cosmological scales. Dark matter is also needed to understand the 
formation of large scale structures in the universe. 

Many candidates have been proposed, either baryonic or not, to explain dark matter. It 
is beyond the scope of this lecture to go through all of these candidates. Here, we restrict 
ourself to the dark matter in galactic halos, and in particular in the halo of our own Galaxy. 
First, we review the evidence for dark matter near the Sun and in the halo of spiral galaxies. 
In Section 3 we discuss the constraint due to the Pauli exclusion principle on neutrinos 
as a candidate for halo dark matter. Based on that argument neutrinos can practically 
be excluded as a major constituent of the dark galactic halos. In Section 4 we present 
the baryonic candidates and in Section 5 we elaborate in some detail on the detection of 
MACHOs (Massive Astrophysical Compact Halo Objects) through microlensing as well as 
on the most recent observations. Section 6 is devoted to a scenario in which part of the dark 
matter is in the form of cold molecular clouds (mainly of H 2 ). 

2 Evidence for dark matter 

In this Section we briefly outline the evidence for dark matter in the solar neighbourhood 
and in the halos of spiral galaxies. Moreover, we discuss also the total mass of our Galaxy, 
which can be inferred from studies of the proper motion of the satellites of the Milky Way. 

2.1 Dark matter near the Sun 

The local mass density || in main sequence and giant stars, stellar remnants (directly ob- 
served or inferred from models of galactic and stellar evolution), gas and dust yields a lower 
limit to the total density of p ~ 0.114 M Q pc~ 3 . Correspondingly, one finds a mass-to-light 
ratio of 

T ~ 1.7T (1) 

(T = Mq/L q , where M® is the mass and L the luminosity of the Sun, respectively). 
The local mass density is determined from carefully selected star samples by analyzing the 
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velocity dispersion and density profile in the direction normal to the galactic plane. This 
yields a total density p = 0.18 ± 0.03 M Q pc~ 3 for the local matter, or equivalently 

T ~ 2.7T . (2) 

Therefore, at least 0.03 M & pc~ 3 is the contribution from dark matter. The presence of disk 
dark matter has long been suspected || and it is most likely baryonic. 

Recently, at least 8 brown dwarfs have been detected within a distance of about 100 
light years from the Sun. One of these brown dwarfs is about 70 million years old and has 
an estimated mass of 0.065 M & . Moreover, some brown dwarfs have been found orbiting 
brighter compagnons, and other as free flying in the Pleiades cluster. It is still premature to 
infer on their contribution to the local dark matter, although it is plausible that they may 
make up an important fraction, if not all. 



2.2 Rotation curves of spiral galaxies 

The best evidence for dark matter in galaxies comes from the rotation curves of spirals. 
Measurements of the rotation velocity v rot of stars up to the visible edge of the spiral galaxies 
and of HI gas in the disk beyond the optical radius (by measuring the Doppler shift in the 
21-cm line) imply that v rot ~ constant out to very large distances, rather than to show 
a Keplerian falloff. These observations started around 1970 ||, thanks to the improved 
sensitivity in both optical and 21-cm bands. By now there are observations for over thousand 
spiral galaxies with reliable rotation curves out to large radii |7J. In almost all of them the 
rotation curve is flat or slowly rising out to the last measured point. Very few galaxies show 
falling rotation curves and those that do either fall less rapidly than Keplerian have nearby 
companions that may perturb the velocity field or have large spheroids that may increase 
the rotation velocity near the centre. 

One of the best exemple for the measurement of the rotation velocity is the spiral galaxy 
NGC 3198 [§] (see Fig. 1). Its halo density can be fitted by the following formula 

p(r) = TTW ' (3) 

where po = 0.013/iq M q pc~ 3 (h Q being the Hubble constant in units of H = 100/io km s _1 
kpc -1 , and 0.4 < ho < 1), a = 6.4/ig 1 kpc, and 7 = 2.1. The total mass inside the last 
measured point of the rotation curve is 1.1 x lO n ho 1 M Q , which yields a total mass-to-light- 
ratio T = 28h T Q . This has to be considered as a lower limit, since there is certainly still a 
lot of dark matter beyond the last measured point on the rotation curve. The dark halo is 
at least four times as massive as the disk. Such a value for the mass-to-light-ratio is typical 
for spiral galaxies. Similar conclusions hold also for elliptical galaxies, although one cannot 
measure rotation velocities. 



2.3 Mass of the Milky Way 

There are also measurements of the rotation velocity for our Galaxy. However, these ob- 
servations turn out to be rather difficult, and the rotation curve has been measured only 
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up to a distance of about 20 kpc. Without any doubt our own galaxy has a typical flat 
rotation curve. A fact this which imply that it is possible to search directly for dark matter 
characteristic of spiral galaxies in our own Milky Way 

In oder to infer the total mass one can also study the proper motion of the Magellanic 
Clouds and of other satellites of our Galaxy. Recent studies || [10], |llj do not yet allow 



an accurate determination of v ro t(LMC)/vo (vq = 210 ± 10 km/s being the local rotational 
velocity). Lin et al. [10| analyzed the proper motion observations and concluded that within 
100 kpc the Galactic halo has a mass ~ 5.5 ± 1 X 10 11 M Q and a substantial fraction ~ 50% 
of this mass is distributed beyond the present distance of the Magellanic Clouds of about 
50 kpc. Beyond 100 kpc the mass may continue to increase to ~ 1O 12 M within its tidal 
radius of about 300 kpc. This value for the total mass of the Galaxy is in agreement with 
the results of Zaritsky et al. |J, who found a total mass in the range 9.3 to 12.5 xl0 11 M Q , 
the former value by assuming radial satellite orbits whereas the latter by assuming isotropic 
satellite orbits. 

The results of Lin et al. [T(J suggest that the mass of the halo dark matter up to the 
Large Magellanic Cloud (LMC) is roughly half of the value one gets for the standard halo 
model (with flat rotation curve up to the LMC and spherical shape), implying thus the 
same reduction for the number of expected microlensing events. Kochanek (Tl| analysed 
the global mass distribution of the Galaxy adopting a Jaffe model, whose parameters are 
determined using the observations on the proper motion of the satellites of the Galaxy, the 
Local Group timing constraint and the ellipticity of the M31 orbit. With these observations 
Kochanek [|11[] concludes that the mass inside 50 kpc is 5.4 ± 1.3 x lO n M . This value 
becomes, however, slightly smaller when using only the satellite observations and the disk 
rotation constraint, in this case the median mass interior to 50 kpc is in the interval 3.3 to 
6.1 (4.2 to 6.8) without (with) Leo I satellite in units of 1O 11 M . The lower bound without 
Leo I is 65% of the mass expected assuming a flat rotation curve up to the LMC. 



3 Neutrinos as halo dark matter 

For stable neutrinos (with mass < 1 MeV) one gets the following cosmological upper bound 
on the sum of their masses |L2], |13| 

m u < 200h 2 eV . (4) 

V 

If neutrinos make up the dark matter in the galactic halos, we may decribe them as forming 
a bound system which resembles in the central regions to an isothermal gas sphere. The core 
radius of such an isothermal sphere is 



1/2 

(5) 



where a is the one-dimensional velocity dispersion and po is the central density. The velocity 
distribution of the neutrinos is Maxwellian and the maximum phase-space density is 

4-5 
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The requirement that the maximum phase-space density does not violate the Pauli exclusion 
principle (n c < g u /h 3 , where g v is the number of helicity states and h Planck's constant) 



leads then to the following lower limit for the neutrino mass [14 



m v > 120 eV 



100 km s^y 74 (i k pc y /2 1/4 



a 



r c J 



(7) 



Typical values for spiral galaxies are o ~ 150 km s _1 and r c ~ 20 kpc. This way we get 
a lower bound m„ > 25 — 30 eV [|TJ|, [15]], which is still consistent with the cosmological 
bound eq.(4). See also refs. ]T7J, 17 1 for a discussion of more precise bounds for spirals and 



ellipticals by considering different visible and dark matter distributions as well as the case 
where the neutrinos are not fully degenerate. However, when considering dwarf galaxies 
for which r c < 2 kpc and a ~ 10 — 30 km s _1 one gets m v > 100 — 500 eV [IS], which 
is clearly in contraddiction with the cosmological bound, excluding thus neutrinos as dark 
matter candidate for the halo of dwarf galaxies and in turn of spiral galaxies. 

This latter point follows also from considering the dwarf galaxies Draco and Ursa Minor, 
which are both satellites of our Galaxy [TOl and, therefore, are located in its halo. In fact, 
if their dark matter halo consist of neutrinos with mass ~ 30 eV, then r c ~ 10 kpc and the 
total mass would be ~ 4 x 10 11 M . However, such a high value for the total mass can 
be excluded by the requirement that the dynamical friction time for such a satellite galaxy 
moving in the halo of our Galaxy has to be longer than the age of Galaxy ~ 10 10 yr. The 
upper value for the total mass one infers this way is of the order of 10 10 M G . Therefore, one 
gets a lower limit of ~ 80 eV for the neutrino mass. 



4 Baryonic dark matter 

Before discussing the baryonic dark matter we would like to mention that another class of 
candidates which is seriously taken into consideration is the so-called cold dark matter, which 



consists for instance of axions or supersymmetric particles like neutralinos pD |. Here, we 
will not discuss cold dark matter in detail. However, recent studies seem to point out that 
there is a discrepancy between the calculated (through N-body simulations) rotation curve 



for dwarf galaxies assuming an halo of cold dark matter and the measured curves |22 
If this fact is confirmed, this would exclude cold dark matter as a major constituent of the 
halo of dwarf galaxies and possibly also of spiral galaxies. 

From the Big Bang nucleosynthesis model |23|, [24| and from the observed abundances 
of primordial elements one infers: 0.010 < IiqQb < 0.016 or with h ~ 0.4 — 1 one gets 
0.01 < Qb < 0.10 (where Qb = PB/Pcrit, and p crit = 3Hq/8ttG). Since for the amount of 
luminous baryons one finds Qi um "C &b, it follows that an important fraction of the baryons 
are dark. In fact the dark baryons may well make up the entire dark halo matter. 

The halo dark matter cannot be in the form of hot ionized hydrogen gas otherwise there 
would be a large X-ray flux, for which there are stringent upper limits. The abundance of 
neutral hydrogen gas is inferred from the 21-cm measurements, which show that its contri- 
bution is small. Another possibility is that the hydrogen gas is in molecular form clumped 



2 One gets a slightly higher bound, by a factor 2 1 / 4 , using the fact that neutrinos behaves practically as 
collisionless particles and thus by applying Liouville's theorem Jli]]. 
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into cold clouds, as we will discuss in some detail in Section 6. Baryons could otherwise 
have been processed in stellar remnants (for a detailed discussion see [25]). If their mass 
is below ~ 0.08 M they are too light to ignite hydrogen burning reactions. The possible 
origin of such brown dwarfs or Jupiter like bodies (called also MACHOs), by fragmentation 
or by some other mechanism, is at present not understood. It has also been pointed out that 
the mass distribution of the MACHOs, normalized to the dark halo mass density, could be 
a smooth continuation of the known initial mass function of ordinary stars [Q. The am- 
bient radiation, or their own body heat, would make sufficiently small objects of H and He 
evaporate rapidly. The condition that the rate of evaporation of such a hydrogenoid sphere 
be insufficient to halve its mass in a billion years leads to the following lower limit on their 
mass [26|: M > IQ~ 7 M & (Ts/30 K) 3 ^ 2 (l g cm~ 3 / p) 1 ^ 2 (Ts being their surface temperature 
and p their average density, which we expect to be of the order ~ 1 g cm~ 3 ). 

Otherwise, MACHOs might be either M-dwarfs or else white dwarfs. As a matter of 
fact, a deeper analysis shows that the M-dwarf option look problematic. The null result of 
several searches for low-mass stars both in the disk and in the halo of our galaxy suggest 
that the halo cannot be mostly in the form of hydrogen burning main sequence M-dwarfs. 
Optical imaging of high-latitude fields taken with the Wide Field Camera of the Hubble 



Space Telescope indicates that less than ~ 6% of the halo can be in this form p7 |. Also 
a substantial component of neutron stars and black holes with mass higher than ~ 1 M 
is excluded, for otherwise they would lead to an overproduction of heavy elements relative 
to the observed abundances. A scenario with white dwarfs as a major constituent of the 



galactic halo dark matter has been explored [p8 ]. However, it requires a rather ad hoc initial 



mass function sharply peaked around 2-6 M . Future Hubble deep field exposures could 
either find the white dwarfs or put constraints on their fraction in the halo |[29||. 



5 Detection of MACHOs through microlensing 



It has been pointed out by Paczyhski |30| that microlensing allows the detection of MACHOs 

7 < M/M G < 



located in the galactic halo in the mass range 
the French collaboration EROS ]31 



10~ Y < M/M < 1. In September 1993 
announced the discovery of 2 microlensing candidates 
and the American-Australian collaboration MACHO of one candidate [32]. In the meantime 
the MACHO team reported the observation of altogether 8 events (one of which is a binary 
lensing event) analyzing 2 years of their data by monitoring about 8.5 million of stars in the 
LMC Their analysis leads to an optical depth of r = 2.91^9 x 10~ 7 an d correspondingly 
to a halo MACHO fraction of the order of 45-50% and an average mass O.5±g;|M , under 
the assumption of a standard spherical halo model. It may well be that there is also a 
contribution of events due to MACHOs located in the LMC itself or in a thich disk of our 

-8 



galaxy, the corresponding optical depth is estimated to be r = 5.4 x 10~ 8 ||33[[ . EROS 
has also searched for very-low mass MACHOs by looking for microlensing events with time 
scales ranging from 30 minutes to 7 days |34|]. The lack of candidates in this range places 
significant constraints on any model for the halo that relies on objects in the range 5 x 10~ 8 < 
M/M & < 5 x 10~ 4 . Similar conclusions have been reached also by the MACHO team fl33 
Moreover, the Polish- American team OGLE |]3"5| , the MACHO [3B| and the French DUO | 
collaborations found altogether more than ~ 100 microlensing events by monitoring stars 
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located in the galactic bulge. The inferred optical depth for the bulge turns out to be higher 
than previously thought. These results are very important in order to study the structure 
of our Galaxy. 

In the following we present the main features of microlensing, in particular its probability 



and the rate of events [3S|. An important issue is the determination from the observations 
of the mass of the MACHOs that acted as gravitational lenses as well as the fraction of halo 
dark matter they make up. The most appropriate way to compute the average mass and 
other important information is to use the method of mass moments developed by De Rujula 
et al. PSA, which will be briefly discussed in Section 5.4. 



5.1 Microlensing probability 

When a MACHO of mass M is sufficiently close to the line of sight between us and a more 
distant star, the light from the source suffers a gravitational deflection (see Fig. 2). The 
deflection angle is usually so small that we do not see two images but rather a magnification 
of the original star brightness. This magnification, at its maximum, is given by 

A U2 + 2 



'-max 



M ( M 2 + 4)1/2 • 



Here u = cI/Re {d is the distance of the MACHO from the line of sight) and the Einstein 
radius Re is defined as 



AGMD 



R E = — x(l-x) (9) 

& 

with x = s/D, and where D and s are the distance between the source, respectively the 
MACHO and the observer. 

An important quantity is the optical depth r opt to gravitational microlensing defined as 

T op t= [ dx^^-p(x)D 2 x(l - x) (10) 
Jo cr 

with p(x) the mass density of microlensing matter at distance s = xD from us along the line 
of sight. The quantity r opt is the probability that a source is found within a radius Re of 
some MACHO and thus has a magnification that is larger than A = 1.34 (d < Re)- 
We calculate r opt for a galactic mass distribution of the form 



Po (a 2 + R? 
a 2 + r 2 



p(ri = ruv i 9 1 z GC > , (n; 



I r | being the distance from the Earth. Here, a is the core radius, po the local dark mass 
density in the solar system and Rgc the distance between the observer and the Galactic 
centre. Standard values for the parameters are po = 0.3 GeV/cm 3 = 7.9 lO _3 M /pc 3 , a = 
5.6 kpc and Rgc = 8.5 kpc. With these values we get, for a spherical halo, r opt = 0.7 x 10~ 6 
for the LMC and r opt = 10" 6 for the SMC g|. 

The magnification of the brightness of a star by a MACHO is a time-dependent effect. For 
a source that can be considered as pointlike (this is the case if the projected star radius at 



6 



the MACHO distance is much less than Re) the light curve as a function of time is obtained 
by inserting 

id 2 + vlt 2 ) 1 / 2 

u(t) = {d +V / } (12) 
Ke 

into eq.([8]), where vt is the transverse velocity of the MACHO, which can be inferred from 
the measured rotation curve (vt ~ 200 km/ s). The achromaticity, symmetry and uniqueness 
of the signal are distinctive features that allow to discriminate a microlensing event from 
background events such as variable stars. 

The behaviour of the magnification with time, A(t), determines two observables namely, 
the magnification at the peak A(0) - denoted by A max - and the width of the signal T (defined 
as being T = R e /vt)- 



5.2 Microlensing rates 

The microlensing rate depends on the mass and velocity distribution of MACHOs. The 
mass density at a distance s = xD from the observer is given by eq. (|TT|) ■ The isothermal 
spherical halo model does not determine the MACHO number density as a function of mass. 
A simplifying assumption is to let the mass distribution be independent of the position in 
the galactic halo, i.e., we assume the following factorized form for the number density per 
unit mass dn/dM, 

dU Alt dn ° a2 + R GC . dn UiT(\j /i q\ 

-dM = — 2 » 2 — — z— dfx = ——H(x)dn , (13) 



dM dfi a 2 + Rq C + D 2 x 2 — 2DR GC xcosa d\i 



with fi = M/Mq (a is the angle of the line of sight with the direction of the galactic centre) 



n not depending on x and is subject to the normalization / rf/i^M = p . Nothing a priori 



is known on the distribution dn^/dM. 

A different situation arises for the velocity distribution in the isothermal spherical halo 
model, its projection in the plane perpendicular to the line of sight leads to the following 
distribution in the transverse velocity 

fM = \ VT e-'»T' v H. (14) 



v 



(vh ~ 210 km/s is the observed velocity dispersion in the halo). 

In order to find the rate at which a single star is microlensed with magnification A > A min , 
we consider MACHOs with masses between M and M + 8M, located at a distance from the 
observer between s and s + 5s and with transverse velocity between vt and vt + 8v T . The 
collision time can be calculated using the well-known fact that the inverse of the collision 
time is the product of the MACHO number density, the microlensing cross-section and the 
velocity. The rate dT, taken also as a differential with respect to the variable u, at which a 
single star is microlensed in the interval djjbdudvrdx is given by |39], [IT 



dT = 2vt f (vt) Dr e[^x(1 — x)] 1 ^ 2 H(x)—^-dfj,dudvTdx, (15) 

(JjUj 
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with 



AGM Q D 



ri = ^ (3.2 x Wkm) 1 . (16) 

cr 

One has to integrate the differential number of microlensing events, dN ev = N^t b s dT, over 
an appropriate range for /i, x, u and vt, in order to obtain the total number of microlensing 
events which can be compared with an experiment monitoring N* stars during an observation 
time t b s and which is able to detect a magnification such that A max > Ath- The limits of 
the u integration are determined by the experimental threshold in magnitude shift, Attith' 
we have < u < uth- 

The range of integration for fi is where the mass distribution dri® / dfi is not vanishing and 
that for x is < x < D^/D where is the extent of the galactic halo along the line of 
sight (in the case of the LMC, the star is inside the galactic halo and thus D^j D = 1.) The 
galactic velocity distribution is cut at the escape velocity v e ~ 640 km/s and therefore vt 
ranges over < vt < v e . In order to simplify the integration we integrate vt over all the 
positive axis, due to the exponential factor in /(vt) the so committed error is negligible. 

However, the integration range of dfidudvxdx does not span all the interval we have just 
described. Indeed, each experiment has time thresholds T min and T max and only detects 
events with: T m i n < T < T max , and thus the integration range has to be such that T lies in 
this interval. The total number of micro-lensing events is then given by 



N, 



cv 



J dN ev Q(T - T mm )Q(T max - T), (17) 



where the integration is over the full range of dfidudvTdx. T is related in a complicated way 
to the integration variables, because of this, no direct analytical integration in eq. (|i~7D can 
be performed. 

To evaluate eq. (|i~7D we define an efficiency function eo(/i) which measures the fraction of 
the total number of microlensing events that meet the condition on T at a fixed MACHO mass 
M = JiM & . A more detailed analysis |39j shows that eo(/^) is in very good approximation 



equal to unity for possible MACHO objects in the mass range of interest here. We now can 
write the total number of events in eq.(|i~7l) as 



N ev = / dN ev e (/x). 



Due to the fact that eo is a function of fi alone, the integration in diidudvxdx factorizes into 
four integrals with independent integration limits. 

In order to quantify the expected number of events it is convenient to take as an example 
a delta function distribution for the mass. The rate of microlensing events with A > A min 
(or u < u max ), is then 

T(A min ) = tu max = Dr E u max ^ vhyj — 7= / dx[x(l - x)] 1/2 H(x) . (19) 

M Q JO 

Inserting the numerical values for the LMC (D=50 kpc and a = 82°) we get 



-13 1 / v h \ ( 1 \ I Po 



T = 4 x 1(T 13 - — - (20) 

a \2W km/s) {^JkicJ \0.3GeV/cm3) ' 1 ' 
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For an experiment monitoring N* stars during an observation time t b s the total number of 



events with a magnification A > A min is: N ev (A % 



1 we show some values of N e 
and A min = 1.34 (or Am min 

Table 1 



for the LMC, taking t obs 
= 0.32). 



N ic t b s T(A min ). In the following Table 
10 7 s (~ 4 Months), = 10 6 stars 



MACHO mass in units of M Q 


Mean Re in km 


Mean microlensing time 


N 

1 > ev 


io- 1 


0.3 x 10 9 


1 month 


1.5 


io- 2 


10 8 


9 days 


5 


io- 4 


10 7 


1 day 


55 


IO" 6 


10 6 


2 hours 


554 



Gravitational microlensing could also be useful for detecting MACHOs in the halo of 



nearby galaxies f42| , [431 such as M31 or M33, for which experiments are under way. In fact, 
it turns out that the massive dark halo of M31 has an optical depth to microlensing which is 
of about the same order of magnitude as that of our own galaxy 



10" 



44 . Moreover, 



an experiment monitoring stars in M31 would be sensitive to both MACHOs in our halo and 
in the one of M31. One can also compute the microlensing rate |44j for MACHOs in the 
halo of M31, for which we get 



1.8 x 10~ 12 - I 

s \210 km/s 




1 



P(0) 



D/kpc) \lGev/cm 3 J Jm/M q 



(21) 



(p(0) is the central density of dark matter.) In the following Table 2 we show some values 
of Nt due to MACHOs in the halo of M31 with t nhq = 10 7 s and N+ = 10 6 stars. In the last 



column we give the corresponding number of events, N ev , due to MACHOs in our own halo. 
The mean microlensing time is about the same for both types of events. 

Table 2 



MHO mass in units of M Q 


Mean Re in km 


Mean microlensing time 


N a 

ev 


N 


io- 1 


7 x 10 8 


38 days 


2 


1 


io- 2 


2 x 10 s 


12 days 


7 


4 


10" 4 


2 x 10 7 


30 hours 


70 


43 


IO" 6 


2 x 10 6 


3 hours 


700 


430 



Of course these numbers should be taken as an estimate, since they depend on the details 
of the model one adopts for the distribution of the dark matter in the halo. 



5.3 Most probable mass for a single event 

The probability P that a microlensing event of duration T and maximum amplification A. 
be produced by a MACHO of mass /i (in units of M ) is given by ||IJ 

P{^T) oc f 4 J l G dx{x{l-x)fH{x)exp (~ ! ^^) , ( 
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which does not dependent on A max and P(/i,T) = P(/z/T 2 ). The measured values for T 
are listed in Table 3, where [Imp is the most probable value. We find that the maximum 
corresponds to [ir 2 E jv 2 H T 2 = 13.0 fl45| , f46[1 . The 50% confidence interval embraces for the 
mass n approximately the range 1/3/zmp up to 3/^mp- Similarly one can compute P(/i, T) 
also for the bulge events (see [fi6j|). 

Table 3: Values of ump (in M Q ) for eight microlensing events detected in the LMC (A4 
= American-Australian collaboration events (i = 1,..,6); Pi and P2 French collaboration 
events). For the LMC: v H = 210 km s" 1 and r E = 3.17 x 10 9 km. 





At 


A 2 


^3 


A, 


A, 


A, 


Pi 


P 2 


T (days) 


17.3 


23 


31 


41 


43.5 


57.5 


27 


30 




0.099 


0.132 


0.177 


0.235 


0.249 


0.329 


0.155 


0.172 




0.13 


0.23 


0.41 


0.72 


0.81 


1.41 


0.31 


0.38 



5.4 Mass moment method 

A more systematic way to extract information on the masses is to use the method of mass 
moments [p9| . The mass moments < fi m > are defined as 

<//">=/ da eM . (23) 

J an 

< /i m > is related to < r n >= J2events rn '■> with r = (t>^/r s )T, as constructed from the 
observations and which can also be computed as follows 

<r n >= J dN ev e n (p) r n = Vu TH T(2 - m)H{m) < fi m > , (24) 

with m = (n + l)/2 and 

V = 2Nj obs Dr E v H = 2.4 x 10 3 pc 3 — f^° bs , (25) 

10° star — years 

T(2 - m) = [°° (^) ^ f(v T )dv T , (26) 
Jo \v H J 

Him) = f (x(l ~ x)) m H(x)dx . (27) 
Jo 

The efficiency e n (/i) is determined as follows (see |39| ) 



_ SdNIM 6(T)r" 
eM = Trf^W r" ' (28) 

where dN* v (p) is defined as dN ev in eq.(^) with the MACHO mass distribution concentrated 
at a fixed mass £1: dn /dfi = no 8{fi — p)/ ' \i. e(T) is the experimental detection efficiency. 
For a more detailed discussion on the efficiency see ref. [[|7j . 
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A mass moment < /i m > is thus related to < r n > as given from the measured values of 
T in a microlensing experiment by 

< r" > 

< /i m >= —-, . (29) 

Vu TH T(2-m)H(m) 

The mean local density of MACHOs (number per cubic parsec) is < fi° >. The average 
local mass density in MACHOs is < fi 1 > solar masses per cubic parsec. In the following we 
consider only 6 (see Table 3) out of the 8 events observed by the MACHO group, in fact the 
two events we neglect are a binary lensing event and an event which is rated as marginal. 
The mean mass, which we get from the six events detected by the MACHO team, is 

1^ = 0.27M S . (30) 

(To obtain this result we used the values of r as reported in Table 3, whereas T(1)H(1) = 
0.0362 and T(2)H(0) = 0.280 as plotted in Fig. 6 of ref. ||). If we include also the two 



EROS events we get a value of 0.26 M for the mean mass. The resulting mass depends on 
the parameters used to describe the standard halo model. In order to check this dependence 
we varied the parameters within their allowed range and found that the average mass changes 
at most by ± 30%, which shows that the result is rather robust. Although the value for the 
average mass we find with the mass moment method is marginally consistent with the result 
of the MACHO team, it definitely favours a lower average MACHO mass. 

One can also consider other models with more general luminous and dark matter distri- 
butions, e.g. ones with a flattened halo or with anisotropy in velocity space [HHf, in which 



case the resulting value for the average mass would decrease significantly. If the above value 
will be confirmed, then MACHOs cannot be brown dwarfs nor ordinary hydrogen burning 
stars, since for the latter there are observational limits from counts of faint red stars. Then 
white dwarfs are the most likely explanation. As mentioned in Section 4 such a scenario has 
been explored recently |28[ . However, it has some problems, since it requires that the initial 
mass function must be sharply peaked around 2 — 6 M & . Given these facts, we feel that 
the brown dwarf option can still provide a sensible explanation of the observed microlensing 
events g|. 

Another important quantity to be determined is the fraction / of the local dark mass 
density (the latter one given by po) detected in the form of MACHOs, which is given by 
/ = Mq/p ~ 126 pc 3 < fi 1 >. Using the values given by the MACHO collaboration 



for their two years data |33| (in particular Uth — 0.661 corresponding to A > 1.75 and an 
effective exposure N*t ty S of ~ 5 x 10 6 star-years for the observed range of the event duration T 
between ~ 20 - 50 days) we find / ~ 0.54, which compares quite well with the corresponding 
value (/ ~ 0.45 based on the six events we consider) calculated by the MACHO group in 
a different way. The value for / is obtained again by assuming a standard spherical halo 
model. 
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Table 4: Values of [imp (in M Q ) as obtained by the corresponding P(fi,T) for eleven 



microlensing events detected by OGLE in the galactic bulge 
te = 1-25 x 10 9 km.) (T is in days as above.) 



30 km s and 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


T 


25.9 


45 


10.7 


14 


12.4 


8.4 


49.5 


18.7 


61.6 


12 


20.9 


T 


0.054 


0.093 


0.022 


0.029 


0.026 


0.017 


0.103 


0.039 


0.128 


0.025 


0.043 


^MP 


0.61 


1.85 


0.105 


0.18 


0.14 


0.065 


2.24 


0.32 


3.48 


0.13 


0.40 



Similarly, one can also get information from the events detected so far towards the galactic 
bulge. The mean MACHO mass, which one gets when considering the first eleven events 
detected by OGLE in the galactic bulge (see Table 4), is ~ 0.29M Q 
discovered during the first year of operation by the MACHO team |31 



From the 40 events 
(we considered only 

the events used by the MACHO team to infer the optical depth without the double lens 
event) we get an average value of O.16M . The lower value inferred from the MACHO 
data is due to the fact that the efficiency for the short duration events (~ some days) is 
substantially higher for the MACHO experiment than for the OGLE one. These values for 
the average mass suggest that the lens are faint disk stars. 

Once several moments < fi m > are known one can get information on the mass distribution 
drio/d/j,. Since at present only few events towards the LMC are at disposal the different 
moments (especially the higher ones) can be determined only approximately. Nevertheless, 
the results obtained so far are already of interest and it is clear that in a few years, due also 
to the new experiments under way (such as EROS II and OGLE II), it will be possible to 
draw more firm conclusions. 



6 Dark clusters of MACHOs and cold molecular clouds 

A major problem which arises is to explain the formation of MACHOs, as well as the nature 
of the remaining amount of dark matter in the galactic halo. We feel it hard to conceive a 
formation mechanism which transforms with 100% efficiency hydrogen and helium gas into 
MACHOs. Therefore, we expect that also cold clouds (mainly of H 2 ) should be present in 
the galactic halo. Recently, we have proposed a scenario |5T| in which dark clusters of 
MACHOs and cold molecular coulds naturally form in the halo at galactocentric distances 
larger than 10-20 kpc, where the relative abundance depends on the distance (similar ideas 



have also been developed in refs. |52|, |53[]). Our scenario can be summarized as follows. 

After its initial collapse, the proto galaxy (PG) is expected to be shock heated to its 
virial temperature ~ 10 6 K. Since overdense regions cool more rapidly than average (by 
hydrogen recombination), proto globular cluster (PGC) clouds form in pressure equilibrium 
with diffuse gas. At this stage, the PGC cloud temperature is ~ 10 4 K, its mass and size 
are ~ 10 6 (R/ 'kpc) 1 / 2 M and ~ 10 (R/kpc) l l 2 pc, respectively. The subsequent evolution of 
the PGC clouds will be different in the inner and outer part of the galaxy, depending on the 
decreasing collision rate and ultraviolet (UV) fluxes as the galactocentric distance increases. 
Below 10 4 K, the main coolants are H2 molecules and any heavy element produced in a 
first chaotic galactic phase. In the central region of the galaxy an Active Galactic Nucleus 
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and/or a first population of massive stars are expected to exist, which act as strong sources 
of UV radiation that dissociates the H 2 molecules present in the inner part of the halo. As a 
consequence, cooling is heavily suppressed and so inner PGC clouds remain for a long time 
at temperature ~ 10 4 K, resulting in the imprinting of a characteristic mass ~ 1O 6 M . Later 
on, the cloud temperature suddenly drops below 10 4 K and the subsequent evolution leads 
to the formation of stars and ultimately to stellar globular clusters. In the outer regions of 
the halo the UV-flux is suppressed, so that no substantial H 2 depletion actually happens. 
This fact has three distinct implications: (i) no imprinting of a characteristic PGC cloud 
mass shows up, (ii) the Jeans mass can now be lower than 10 _1 M Q , (m) the cooling time is 
much shorter than the collision time. PGC clouds subsequently fragment into smaller clouds 
that remain optically thin until the minimum value of the Jeans mass is attained, thus 
leading to MACHO formation in dark clusters. Moreover, because the conversion efficiency 
of the constituent gas in MACHOs could scarcely have been 100%, we expect the remaining 
fraction of the gas to form self-gravitating molecular clouds, since, in the absence of strong 
stellar winds, the surviving gas remains bound in the dark cluster, but not in diffuse form 
as in this case the gas would be observable in the radio band. 



6.1 Observational Tests 

Let us now address the possible signatures of the above scenario, in addition to the single 
MACHO detection via microlensing. 

We proceed to estimate the 7-ray flux produced in molecular clouds through the interac- 
tion with high-energy cosmic-ray protons. Cosmic rays scatter on protons in the molecules 
producing 7r°'s, which subsequently decay into 7's. An essential ingredient is the knowledge 
of the cosmic ray flux in the halo. Unfortunately, this quantity is experimentally unknown 
and the only available information comes from theoretical estimates. More precisely, from the 
mass-loss rate of a typical galaxy we infer a total cosmic ray flux in the halo F ~ 1.1 x 10~ 4 
erg cm -2 s" 1 . We also need the energy distribution of the cosmic rays, for which we assume 
the same energy dependence as measured on the Earth. We then scale the overall density in 
such a way that the integrated energy flux agrees with the above value. Moreover, we assume 
that the cosmic ray density scales as R~ 2 for large galactocentric distance R. Accordingly, 



we obtain p0|, 51 



$ CR (E, R) ~ 1.9 x 10- 3 ** R (E) a2 a2 + + R £c , (31) 

where <&® R (E) is the measured primary cosmic ray flux on the Earth, a ~ 5 kpc is the halo 
core radius and Rqc ~ 8.5 kpc is our distance from the galactic center. The source function 
g 7 (r), which gives the photon number density at distance r from the Earth, is 

47r r 

? 7 (r) = —pH 2 {r) J dE p $ C R(E p ,R(r)) a iri {p lab ) < n^(E p ) > . (32) 

Actually, the cosmic ray protons in the halo which originate from the galactic disk are mainly 
directed outwards. This circumstance implies that the induced photons will predominantly 
leave the galaxy. However, the presence of magnetic fields in the halo might give rise to 
a temporary confinement of the cosmic ray protons similarly to what happens in the disk. 
In addition, there could also be sources of cosmic ray protons located in the halo itself, 
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as for instance isolated or binary pulsars in globular clusters. As we are unable to give 
a quantitative estimate of the above effects, we take them into account by introducing an 
efficiency factor e, which could be rather small. In this way, the 7-ray photon flux reaching 
the Earth is obtained by multiplying g 7 (r) by e/4irr 2 and integrating the resulting quantity 
over the cloud volume along the line of sight. 

The best chance to detect the 7-rays in question is provided by observations at high 
galactic latitude. Therefore we find 

$ 7 (90°) ^ ef 3.5 x ID- . (33) 

cm 2 s sr 

The inferred upper bound for 7-rays in the 0.8 - 6 GeV range at high galactic latitude 
is 3 x 10~ 7 photons cm~ 2 s _1 sr -1 j54]. Hence, we see from eq. (^) that the presence of 
halo molecular clouds does not lead nowadays to any contradiction with such an upper limit, 
provided ef < 1CT 1 . 

Molecular clouds can be detected via the anisotropy they would introduce in the Cosmic 
Background Radiation (CBR), even if the ratio of the temperature excess of the clouds to 
the CBR temperature is less than ~ 10~ 3 . Consider molecular clouds in M31. Because we 
expect they have typical rotational speeds of 50 — 100 km s -1 , the Doppler shift effect will 
show up as an anisotropy in the CBR. The corresponding anisotropy is then p5f 



^ = ±- S f r v , (34) 
l r c 

where S is the spatial filling factor and T r is the CBR temperature. If the clouds are optically 
thick only at some frequencies, one can use the average optical depth over the frequency range 
of the detector f . We estimate the expected CBR anisotropy between two fields of view (on 
opposite sides of M31) separated by ~ 4° and with angular resolution of ~ 1°. Supposing 
that the halo of M31 consists of ~ 10 6 dark clusters and that all of them lie between 25 kpc 
and 35 kpc, we would be able to detect 10 3 — 10 4 dark clusters per degree square. Scanning 
an annulus of 1° width and internal angular diameter 4°, centered at M31, in 180 steps of 1°, 
we would find anisotropies of ~ 2 x 10 -5 / f in AT/T r (as now S = 1/25). In conclusion, 
the theory does not permit to establish whether the expected anisotropy lies above or below 
current detectability (~ 10~ 6 ), and so only observations can resolve this issue. 

An attractive strategy to discover the halo molecular clouds clumped into dark clusters 



relies upon the absorption lines they would introduce in the spectrum of a LMC star |hj • 

Let us now turn to the possibility of detecting MACHOs in M31 via their infrared emission. 
For simplicity, we assume all MACHOs have equal mass ~ 0.08 M (which is the upper mass 
limit for brown dwarfs) and make up the fraction / of the dark matter in M31. In addition, 
we suppose that all MACHOs have the same age t ~ 10 10 yr | 57[| . As a consequence, 



MACHOs emit most of their radiation at the wavelength X ma x ~ 2.6 fxm. The infrared 
surface brightness I v (b) of the M31 dark halo as a function of the projected separation b 
(impact parameter) is given by 



x 3 a 2 f L 2 — b 2 

iJb) ~ 5 x 10 5 arctanW— — — Jy sr -1 , (35) 

V ; e* - 1 DVa 2 + b 2 V a? + b 2 3 ' K ' 
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where the M31 dark halo radius is taken to be L ~ 50 kpc. Some numerical values of Iv max {b) 
with b = 20 and 40 kpc are ~ 1.6 x 10 3 / Jy sr _1 and ~ 0.4 x 10 3 / Jy sr" 1 , respectively. 
The planned SIRTF Satellite contains an array camera with expected sensitivity of ~ 1.7 x 
10 3 Jy sr _1 per spatial resolution element in the wavelength range 2-6 /im. Therefore, the 
MACHOs in the halo of M31 can, hopefully, be detected in the near future. 

7 Conclusions 

The mistery of the dark matter is still unsolved, however, thanks to the ongoing microlensing 
experiments there is hope that progress on its nature in the galactic halo can be achieved 
within the next few years. It is well plausible that only a fraction of the halo dark matter 
is in form of MACHOs, either brown dwarfs or white dwarfs, in which case there is the 
problem of explaining the nature of the remaining dark matter and the formation of the 
MACHOs. Before invoking the need for new particles as galactic dark matter candidates for 
the remaining fraction, one should seriously consider the possibility that it is in the form 
of cold molecular clouds. A scenario this, for which several observational tests have been 
proposed, thanks to which it should be feasible in the near future to either detect or to put 
stringent limits on these clouds. 

I would like to thank F. De Paolis for carefully reading the manuscript. 
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Figure Captions 



Fig. 1: Rotation curve for NGC 3198 according to van Albada et al. ||. The dotted line 
with error bars refers to the optical and 21 cm hydrogen data, while the solid lines 
are theoretical fits. 

Fig. 2: Definition of various quantities describing a microlensing event. The observer is 
O, the source is S and M is the MACHO. 
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